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Abstract. The hypermagnetic helicity density at the electroweak phase transition (EWPT) 
exceeds many orders of magnitude the galactic magnetic helicity density. Together with pre- 
vious magnetic helicity evolution calculations after the EWPT and hypermagnetic helicity 
conversion to the magnetic one at the EWPT, the present calculation completes the descrip- 
tion of the evolution of this important topological feature of cosmological magnetic fields. It 
suggests that if the magnetic field seeding the galactic dynamo has a primordial origin, it 
should be substantially helical. This should be taken into account in scenarios of galactic 
magnetic field evolution with a cosmological seed. 
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1 Introduction 

The magnetic helicity is a relatively new and very attractive point of interest in cosmic 
magnetohydro dynamics (MHD) and dynamo theory. The point is that the magnetic helicity 
H = f ABd 3 x where B is the magnetic field and A is the vector potential is an integral 
of motion in MHD in the absence of viscosity (inviscid case). It looks natural to base an 
understanding of cosmic MHD on balance equations for conserved quantities such as magnetic 
field energy and magnetic helicity. It has been noted that the magnetic helicity conservation 
is much more restrictive in astrophysical objects than the energy conservation [1, 2]. The 
point is that usually there is a huge supply of kinetic energy in the form of a general rotation 
of a celestial body and it is quite easy to imagine a spectral energy flux which influences 
energy balance including magnetic field energy. 

The nature of the initial fields (and corresponding magnetic helicities) that seed sub- 
sequent dynamo or turbulent amplifications is largely unknown [3, 4]. It might be that the 
seed fields are produced during epoch of galaxy formation from frozen-in magnetic fields of 
protogalaxy experiencing gravitational collapse, or ejected by the first supernovae or active 
galactic nuclei. Let us call this as the standard astrophysical scenario A. Alternatively the 
seed fields might originate from much earlier epochs of the Universe expansion, down to the 
cosmological inflation phase transition epoch [5]. Let us call that as the cosmological scenario 
B. The standard view in scenario A is that the magnetic field evolution in the early Universe 
(or an astrophysical object) starts from a state with almost vanishing magnetic helicity. The 
dynamo amplification of a seed field in astrophysical scenarios produces, however, a large- 
scale magnetic field with substantial magnetic helicity. In such case one must compensate it 
providing a contribution of magnetic helicity of small-scale magnetic fields. Then one faces 
with a severe problem of how to redistribute magnetic helicity over the desired scales in or- 
der to keep it small. Moreover, it is difficult to prevent a disastrous dynamo suppression of 
large-scale magnetic helicity by helical small-scale magnetic field. Alternatively in scenario 
B, if the seed magnetic field contained a lot of magnetic helicity one can use that as the he- 
licity required to obtain the desired large-scale field at later epochs. This makes the dynamo 
generation of galactic magnetic fields much less constrained than in the standard scenario A. 
If the galactic dynamo exploits the primordial magnetic helicity then one must expect that 
a large-scale galactic magnetic field has also a preferable sign of helicity. Observations show 
indeed some hint that one of the possible helicity signs seems preferable [6]. 
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Relying on scenario B with the hypermagnetic field evolution passing through the elec- 
troweak phase transition (EWPT) epoch, we explore here the magnetic helicity generation 
in the early universe. This suggests as a new alternative the possibility that the magnetic 
field starts from a state with a substantial supply of primordial magnetic helicity. Indeed, a 
very small fluctuation of a seed hypermagnetic field at very early epochs before the EWPT 
woul be exponentially amplified due to the presence of a large right electron asymmetry, 
£ e R = fien/T / [7-9] and hence acquire a huge initial magnetic helicity by the EWPT 
epoch and subsequently. This is possible only for hypermagnetic fields having a non-trivial 
topological structure with non- vanishing linkage number n / [10]. The hypermagnetic 
helicity density given by the product hy = Y ■ By, hence proportional to a large hypermag- 
netic field squared value hy ~ By A- This gets transformed to the magnetic helicity, which 
can therefore be much larger than the helicity density associated to galactic magnetic fields, 
hy » /igai where h ga \ ~ £g a ^gal- For £? ga i ~ 10~ 6 G, Z ga i = 100 kpc = 3 x 10 23 cm, this is 
estimated as h gai \ = 3 x 10 11 G 2 cm. Thus, the primordial magnetic helicity can be the main 
supply for magnetic helicity of galaxies. The goal of this paper is to provide a complete de- 
scription of magnetic helicity evolution passing through the EWPT epoch (we also comment 
on later epochs within the causal picture). In Section 2 we calculate the hypermagnetic he- 
licity in the symmetric phase of the hot primordial plasma, using the corresponding solution 
of the Faraday equation. In Section 3 we review the evolution of the hypermagnetic helicity 
passing through the EWPT and in Section 4 we make the final comments on our results and 
on the perspectives for the cosmological origin of helical galactic magnetic fields. 



2 Hypermagnetic helicity 

In the comoving frame V = 1 the Faraday induction equation governing the evolution of 
hypermagnetic fields By = V x Y reads 

(9B 

-^=Vx«yBy+7 |F V%, (2.1) 

where the hypermagnetic helicity coefficient ay [9] is given by the right electron chemical 
potential fx e u and the hot plasma conductivity cr con d{T) ~ T as 

and rjy = (ccond) -1 is the hypermagnetic diffusion coefficient, g = e/ cos^ is the Standard 
Model U(l) gauge coupling. We assume that a certain right electron asymmetry ~ ^ e R(to) at 
a very early cosmological epoch , to £ew, has been generated by an unspecified mechanism. 
This is our starting point. Then in the presence of the hypercharge field Yn this asymmetry 
n eR — n eR = HeRT 2 evolves due to the Abelian anomaly for right electrons, 

d,f eR = - 9 ^y^, y R = -2, (2.3) 

which evolves together with the hypermagnetic field in Eq.(2.1) in a self-consistent way. Note 
that such coupled evolution of By(t) and /i e i?(i) has been recently considered in Ref. [11] for 



1 Note that common expansion can be easily taken into account via conformal coordinates with the change 
of the cosmological time t to the conformal one dt —¥ a{t)dr). 
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the particular case of the Chern-Simons wave hypermagnetic field configuration, but without 
considering such important feature as the hypermagnetic helicity which we discuss here. 

Our second assumption is the presence of a non-zero initial hypermagnetic field Bq ^ 0. 
This should be a mean field with a small amplitude provided by some stochastic distribution 
of hypermagnetic fields. 

The key parameter in the master Eq. (2.1) is the helicity parameter for the hyper- 
magnetic field, given in eq. (2.2). This can be obtained from the Chern-Simons term in the 
effective Standard Model Lagrangian density for the hypercharged field : 

Cos = ~^^B Y Y. (2.4) 

An recent interpretation of the Chern-Simons anomaly parameter a Y as a polarization effect 
has been given in Ref. [9] using standard statistical averaging of the right electron pseudovec- 
tor current < e^jj^e > in the Standard Model Lagrangian in vacuum (alternative one-loop 
level calculations in finite temperature field theory were used in [12, 13]). 

Multiplying Eq. (2.1) by the corresponding vector potential and adding the analogous 
expression obtained from the evolution equation governing the vector potential (multiplied 
by hypermagnetic field) and integrating over space, one gets the evolution equation for the 
hypermagnetic helicity Hy = J d 3 xY • By as 

dHy 



dt 



-2 j (Ey ■B Y )d 3 x- j)[Y B Y + 
+Ey x Y]d 2 S = -2r] Y (t) J d 3 x(V x By) • By + 
+2a Y (t) / d 3 xB Y (t). (2.5) 



Note that we have omitted in the last equality the surface integral §(...) since fields vanish at 
infinity during the symmetric phase. However, such surface integral can be important at the 
boundaries of different phases at the electroweak phase transition, T ~ Tew- I n Ref- [10] the 
authors have studied how the hypermagnetic helicity flux penetrates the surface separating 
the symmetric and broken phases, and how the hypermagnetic helicity density h Y = ByY 
converts into the magnetic helicity density h = BA at the EWPT time, see also discussion 
in Sec. 3 below. 

Notice also that the evolution equation in eq. (2.5) is similar to eq. (7) in Ref. [14], 
which holds after the electroweak phase transition, T <C Tew- There the point- like short- 
range Fermi neutrino-plasma interaction mediated by heavy W, Z-bosons was used, instead of 
the long-range interaction through the massless hypercharge field in the unbroken phase. 

Now using the simplest solution of the Faraday equation (2.1) in the a 2 -dynamo [15] 
that corresponds to maximum hypermagnetic field amplification rate 2 , 



B Y {t) = Bq exp 



4r?y(0 



(2.6) 



"This is the case for the particular hypermagnetic field scale A = k = kt)y /ay where k, — 2, k being the 
Fourier wave number in By(x,<) = J(d 3 fe/(27r) 3 )B y (k,t)e ikx . 
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we obtain from Eq. (2.5) the hypermagnetic helicity 



H Y = H Y (t ) + 2(B%) 2 
a Y (t 



d a x / dt a Y {t ) x 
to 



x exp 



to 



Ht") 



dt 



(2.7) 



where Hy(to) is the initial helicity value at the moment to an d we have omitted the hy- 
permagnetic diffusion term. The latter is usually neglected for an ideal Maxwellian plasma 
<7cond — > oo, tjy — > 0, while here we must compare in Eq. (2.5) the second helicity generation 
term ay ~ (ccond) -1 given by Eq. (2.2) with the first diffusion term r\y ~ (<T C ond) _1 - Notice 
that in Ref. [14] one could neglect the diffusion term for the Maxwellian plasma in the broken 
phase, T <C Tew-, since for that case the magnetic helicity coefficent a [14, 16] did not depend 
on the conductivity. 

The direct estimate of the relative magnitude of the terms in the r.h.s. of Eq. (2.5) 
allows us to neglect the diffusion term only for hypermagnetic field inhomogeneity scales 
obeying the inequality: 

(2.8) 



A > — . 



ay 



For an arbitrary (large) scale in our causal scenario, Ijj > A = kt)y jay 3> rjy/ay, the 
amplification of hypermagnetic field is given by [9] 



B Y (t) = B% exp 



Bq exp 



S3 







[(=- 






r xo dx' 




J X x 


given 


by 




(2.9) 



and that of the hypermagnetic helicity is given by 

Hy(t)=Hy(t ) + 2(B^) 2 J d 3 X 

a Y ((t 



dt ay(t ) x 



x exp 





-) 




[(=- 


K? ) 





H Y (t ) + 2{B£ Y I d d x I dtay(t)x 




x exp 



, 1 1 

166 | ? 



"2 



0.0001 



(2.10) 



Here the ratio x = T/Tew = {tEw/t) 1 ^ 2 is given by the Friedman law and = ^ e R/T is 
the dimensionless right electron asymmetry. Thus, from Eq. (2.9) for the extremum value 
K = 2 we obtain the strongest amplification eq. (2.6) and from Eq. (2.10) one finds the 
corresponding value of the hypermagnetic helicity in eq. (2.7). 
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Let us comment on our reference value choice £, e R ~ 10~ 4 used in eqs. (2. 9), (2. 10). 
Taking into account the right electron (positron) asymmetry n e R — ngR = [i e nT 2 /6 in the 
presence of chirality flip processes with the rate T, and substituting the hyperelectric field 
Ey = —V x By + t)y V x By — QyBy [9] into the Abelian anomaly Eq. (2.3) rewritten in 
uniform medium as dt(n e R — ngR) = —{g 2 /4-/r 2 )EyBy one finds the kinetic equation for pL e R: 

dfi eR 6</ 2 (V X By) • By 

~dT = 4vr 2 rv cond " {Tb + r W (2 - n) 

Here the rate Tb = 6(5 2 j Att 2 ) 2 By /T 2 a C ond coming from the helicity term ~ ay occurs in 
strong hypermagnetic fields much bigger than the chirality flip rate r, Tb S> T [9]. 

Under the assumption of slowly changing hypermagnetic fields By(t) ~ const, and 
choosing the Chern-Simons wave configuration of the hypercharge field as, 

Y = Y z = 0, Y x = Y(t)smk z, Y y = Y(t) cos k z, 

for which (V x By) - By = B Y (t)ko, By(t) = koY(t), we can easily solve the kinetic equation 
(2.11) getting: 



£e.R(V 



r B + r 



e -(r B +r)(t-to) + 



Q „ Q „ Q r919 , 
+ r B + r~ r B + r~ r B ~ t 9 ' 2 ' [ZA2) 

Here we used notation Q = -(6g 2 /4ir 2 T 3 cj cond )(V xB y )-B y = -{<6g 2 /^ 2 T^a cond )B 2 k . In 
obtaining (2.12) we neglected the time dependence for times t — > t^w f° r which TstEW ^ 1- 
On the other hand, retaining the time term for the zero initial asymmetry £, e R(to) = we 
get from (2.12) the asymptotical growth of | £ e R \ in a strong hypermagnetic field due to the 
Abelian anomaly: 

1 - o " (r 



r B + r 

~ Tg' 2 I 1 \ ■ 

Taking into account for the survival condition of the Chern-Simons wave versus ohmic 
diffusion, ko < 10~ 7 T, substituting weak coupling g 2 = 0.12 we get the estimate of the 
lepton asymmetry in a strong hypermagnetic field, | £ e R |~ 3 x 10~ 5 , hence we adopted 
£ e /j ~ 10 -4 as the reference value in Eqs. (2. 9), (2. 10) above. For a topologically non-trivial 
3D-hypermagnetic field configuration with linkage (Gauss) number n S> 1 for which the 
pseudoscalar By • (V x By) ~ n one can expect the right electron asymmetry at the level 
^ e j? ~ 10 -4 , which will be used below as an estimated reference value with respect to which 
we choose to normalize the right electron asymmetry £ e #. 

3 Hypermagnetic helicity evolution 

Let us now turn to the evolution of the hypermagnetic (magnetic) helicity through various 
stages in the evolution of the universe, as illustrated in Fig. 1. 
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3.1 Hypermagnetic helicity evolution through the electroweak phase transition 

Let us note that for Higgs masses mu > 80 GeV the electroweak phase transition cannot be 
first order in the minimal standard electroweak theory, so that a smooth cross-over between 
symmetric and broken phases is more likely [17], given the experimental lower bound on 
Higgs masses and the recent hints from the LHC [18]. 

However, in the presence of strong hypermagnetic fields in the primordial plasma the 
dynamics of the phase transition changes in analogy with the superconductivity in the pres- 
ence of magnetic fields: the second order phase transition may become first order [7]. We 
rely here on such scenario assuming the presence of strong hypermagnetic fields for which 
first order EWPT becomes allowed in the mass region 80 GeV < mu < 160 GeV (see Fig. 
8 in [7]) indicated by current LHC data. 

We first consider what happens with the hypermagnetic helicity passing-through the 
electroweak phase transition. For that let us mention results from paper [10] where the 
flow of the hypermagnetic helicity in the embryo of the new (broken) phase was considered. 
If a single bubble of broken phase appears at T = Tew growing with constant velocity, 
R(t) = v(t — tEw) 3 i then a value of the flow of hypermagnetic helicity density through the 
bubble surface is determined by the surface integral in Eq. (2.5) which we neglected above 
for the symmetric phase with boundary at infinity. The result (Eq. (17) in [10]) shows that 
the value of the hypermagnetic helicity density penetrating the surface of a single bubble, 

h Y {t) _ 5xl0~ 3 n ( B Y (t E w) \ 2 ( t-t EW \ 2 , , 

G 2 cm d(cm) \ 1 G J \ t EW J ' 1 ' ' 

is also large. Here the integer n = =pl,=p2,... denotes the number of pairs of linked hyper- 
magnetic field loops (or knot number) for the non-trivial 3D-configuration, as in Eq. (7) in 
Ref. [10]. Note that n is the pseudoscalar entering the Gauss integral for magnetic helicity 

H(t) = J d 3 xh(x,t) = n$i$2 

which changes the sign, n — > —n, after one of the overlapping oriented loops in a pair of 
magnetic closed tubes changes the direction. 

In order to avoid screening of the hyperelectric field Ey and the time component Yo 
over the surface of the bubble the thickness d of the domain wall separating the two phases 
should be less than the Debye radius, d <r£> = y^T^^y /47re 2 n e ~ 10 /Tew : which allows us 
to estimate the factor d~ l in eq. (3.1) as [d(cm)]^ 1 > 10 15 /2. Then substituting the value 
of the hypermagnetic field By(tEw) estimated in in the leptogenesis scenario as By{tEw) ~ 
5 x 10 17 G [9, 21] , one gets from Eq. (3.1) the helicity density 

h/G 2 cm = 6.25 x 10 47 n[(t - %)/%] 2 . 

Such huge value is estimated at the moment of the growth of a bubble of the new phase, e.g. 
for R(t)/lu < [(t — tEw/tEw)] ~ 10~ 7 4 . Taking into account the subsequent conservation 
of the net global helicity, summed over different protogalactic scales, one finds values that 
are much larger than the helicity density associated to galactic magnetic field strengths, 

h g ai ~ 10 11 G 2 cm. 
3 Here time is fixed at Tew, (t — tEw) /tEW < 1, D = 0.1 — 1 accordingly [19, 20]. 

4 Such bubble size is relevant before percolation (collision and following junction) of the two bubbles, see 
Eq. (21) in Ref. [20]. 
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During the electroweak phase transition the hypermagnetic helicity density converted 
from the symmetric phase to the Maxwellian one is redistributed within a bubble over its 
volume in correspondence with the ratio of volumes 

Vsurface/Vball = 3d/R(t). 

This is because we calculated the helicity flux density in eq. (3.1) only within a thin spherical 
layer V sur f ace = A7rR 2 (t)d. Then assuming that all bubbles are tight-fitting each other within 
the horizon size Ih^ew) = 1-44 cm or the mean magnetic helicity density coincides with 
that within one bubble, and using the ratio 3d/R(t) we obtain from (3.1) at T ~ Tew the 
magnetic helicity density in the new (broken) phase 



h(x ~ 1) 1.5 x W~ 2 n f B Y (t E w) \ 2 
G 2 



x 



( R(t) \ _ 1.5 x 10~ 9 n (B%^ 2 



\ln(t E w)J 1-44 \IG 

2 



X 



x exp 



166(k - 1) f°° dx (i eR {x) 



x x 2 vo.oooi 



(3.2) 



where we have substituted R/Ih = 10 7 as an estimate of the beginning of percolation 
(junction) of bubbles (see Eq. (21) in Ref. [20]). 

We now turn to a brief discussion of the bounds on the topological linkage number 
n = ±1, ±2, ... (|n| > 1) in eq. (3.2). 

Starting from Gauss formula, 

H = n<& 2 = n5 VA 4 , 

and substituting the corresponding helicity density h = 3H/4itR s , and using the maximum 
helicity density /i max = B 2 A one can find a bound on "n" from the requirement that h < h max : 

4 //A 3 4 x 10 6 { Z eR (TEw) \ 3 

n< 3^UJ = ^^{-^r) ' (3 - 3) 

where we have substituted R = 10~ 7 Ih(Tew) and A = kt]y jcty = 3.3x W 6 k/[Tew (Cefl/0-0001)]. 

On the other hand, from the same bound h < h max using Eq. (3.2) and cancelling 
B 2 (tEw) one finds: 

h(x = 1) = 1.5 x 10~ 9 n(£ e fl(?W)/0.0001) 
B 2 A 2.88 x 3.3 x 10 6 k x 10" 16 < ' 

or 

k > 1.6n ] . (3.4) 

V o.oooi / v ' 

Combining (3.3) and (3.4) from the chain of inequalities we get, 

^ 4x 10 6 fZe R (T E w)\ 3 ^ 4x 10 6 
n < 3 — — < 



0.0001 J (1.6) 3 n 3 ' 
hence we find an upper bound on the linkage number 

| n |< 33 . (3.5) 
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Note that it does not depend on a right electron chemical potential nor on the bubble size 
before percolation, chosen in Eq. (3.3) as R = W~ 7 Ih(Tew) [19]. Indeed, the scale A ~ k is 
proportional to the bubble size R, A ~ R, as seen from Eq. (3.4) where k ~ R through the 
helicity h(x = 1) in Eq. (3.2). 

3.2 Hypermagnetic helicity evolution before the electroweak phase transition 

For an arbitrary scale A = Kr/y/ay, k > 1, we obtain from eq. (2.10) the hypermagnetic 
helicity in the unbroken phase (x > 1) as 

hy(x) _ 
G 2 cm 



-0.88 x 10" 



x cxp 



166(k 



1 G 
1 



dx I £ eR (x ) 
x' 3 \ 0.0001 



dx / £ e i?(x' 



0.0001 



(3.6) 



where we put xq = oo, substituted ay(t) from eq. (2.2), and used the expansion time t = 
(Mpi/l,m^Jlf)/2T 2 = M /2T 2 , omitting the initial helicity value Hy(t ). 

Dividing the helicity density (3.6) by its value at Tew (3-2) we get the ratio valid at 
x > 1: 



hy(x) 
h(x = 1) 



8L5 
n 



dx ( ^eR(x') 



x exp 



166(k - 1) 



0.0001 

2 



i eR {x )\ dx' 



0.0001 



(3.7) 



that for a constant value of the right electron asymmetry £ e #/0.0001 = (3 = const equals to 



hy(x) 8.5 | p 



h(x = 1) 



1 

^2 



exp[-a/3 2 (l - l/x)} 



exp[-a/3 2 (l - l/x)} - exp(-a/3 2 
^2 



(3. 



where a = 166(k — 1)/k 2 and we took into account the negative sign of the right electron 
asymmetry, j3 < 0, as seen, e.g., from Eq. (2.12). For a small parameter aj3 2 <C 1, or 
k S> 166/3 2 that is allowed for larger bubbles R S> 10 7 Z// one gets from (3.8) (using relation 



t/tEW = (Tew/TY 



') the temporal dependence 



ew , 



n 



hy{t < t EW ) _ 4.25 | p 
h(t 

3 



(1 - a/3 2 ) 



tEW 



\3/2 

J +0((c 



tEW 



(3.9) 



which is almost linear in the region t < t^w- Thus, we see from (3.9) that like for EWPT of 
the first order there is the jump of helicity density (here for large scales k 3> 166/3 2 ): 

hy(t < t EW ) 4.25 | P | 



h(tEw) 



n 



(3.10) 
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where the linkage number | n |> 1 has the upper bound given by Eq. (3.5). Note that 
the straight lines in Fig. 1 shown for some particular parameters: a given by k = 2, 10 3 ; 
| j3 |= 0.1, 1; n = 1, 30 correspond at t < tEW to the limiting case presented in Eq. (3.9). 

Let us stress that our assumption of a constant right electron asymmetry £ e #/0.0001 = 
(3 = const is provided by the adiabatic condition dt£, e R ~ which we tacitly assume. 

4 Discussion 

The lepton asymmetry plays a crucial role in the amplification of hypermagnetic fields, which 
become helical and supply magnetic helicity to the cosmological and subsequently to the 
galactic magnetic fields. The lepton asymmetry itself evolves due to the Abelian anomaly for 
the right electron current given by Eq. (2.3) in the unbroken phase of the Standard Model. It 
starts from an initial £, e R,(to) generated by an unspecified mechanism and then rises driven by 
ay -helicity parameter (2.2) arising from the Chern-Simons term (2.4). One should note that 
the absence of the Chern-Simons term in the broken phase [12] does not mean that parity 
violation in electroweak interactions disappears. Indeed the polarization effect leading to the 
helicity a parameter governing Maxwellian field evolution exists due to paramagnetism of 
fermions populating the main Landau level [16]. 

The magnetic helicity parameter a which governs the evolution of the magnetic helicity 
after the electroweak phase transition is shown in Fig. 1 by short horizontal lines for t > tEW- 
This change of the magnetic helicity density profile h(t) is explained by a negligible value 
of the helicity parameter a for weak interactions in broken phase at T <C Tew> 
[14, 16], comparing with ay given by Eq. (2.2) for symmetric phase, a <C ay . 

The jump of helicity density at t = tEW is the topological effect of a difference between 
the volume helicity density entering the first line in Eq. (2.5) and the surface helicity term in 
the same equation. While the volume term gives the smooth function (3.6) in the numerator 
of the ratio in Eq. (3.7) the surface term leads to the helicity density at EWPT ~ n given 
by (3.2) in the denominator of Eq. (3.7). 

Note also that, following standard practice, we have neglected turbulence due to a non- 
zero plasma vorticity arising e. g. through bubble collisions during the electroweak phase 
transition. Such simplification is justified in the treatment of hypermagnetic helicity since 
the fluid velocity V(x, t) does not contribute to helicity evolution. As a result here we have 
confined our attention only to the a 2 -dynamo mechanism, avoiding aVt -dynamo scenario. 

In summary, we have found that the magnetic helicity which becomes an inviscid in- 
variant after the EWPT, varies dramatically before this phase transition. Due to this, the 
cosmological magnetic field becomes helical before the phase transition and remains helical 
after it, at least within the range of applicability of our causal scenario. It means that the 
seed magnetic field for the galactic dynamo if provided by a primordial cosmological magnetic 
field should be substantially helical. This seed magnetic helicity must be taken into account 
in scenarios of galactic magnetic field evolution with a cosmological seed. In particular, the 
intergalactic magnetic field suggested in [22] is expected to be substantially helical. 
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the AHEP group of IFIC for hospitality. This work was supported by the Spanish MEC 
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Figure 1. Magnetic helicity evolution in the EWPT epoch for various parameter choices: top: 
a = 0.166, | (3 |= 1, n = 1; middle: a = 41.375 (i.e. k = 2), | ft |= 0.1, n = 30; bottom: red line 
a = 41.375 , | /3 |= 0.1, n = 1, blue a = 41.375, | /3 \= 0.1, n = 30, green: a = 0.166 (i.e. k = 10 3 ), 
| /3 |= 0.1, n = 1; for t > tsw red and green lines are identical. 
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